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We numerically study turbulent Taylor-Couette flow (TCF) between two independently 
rotating cylinders and the transition to rotating plane Couette flow (RPCF) in the limit 
of infinite radii. By using the shear Reynolds number Res and rotation number as 
dimensionless parameters, the transition from TCF to RPCF can be studied continuously 
without singularities. Already for radius ratios t] ^ 0.9 we find that the simulation 
results for various radius ratios and for RPCF collapse as a function of Rq, indicating a 
turbulent behaviour common to both systems. We observe this agreement in the torque, 
mean momentum transport, mean profiles, and turbulent fluctuations. Moreover, the 
central profiles in TCF and RPCF for Rq > 0 are found to conform with inviscid neutral 
stability. Intermittent bursts that have been observed in the outer boundary layer and 
have been linked to the formation of a torque maximum for counter-rotation are shown to 
disappear as 77 —^ 1. The corresponding torque maximum disappears as well. Instead two 
new maxima of different origin appear for ry ^ 0.9 and RPCF, a broad and a narrow one, 
in contrast to the results for smaller ry. The broad maximum at Rq = 0.2 is connected with 
a strong vortical flow and can be reproduced by streamwise invariant simulations. The 
narrow maximum at Rq — 0.02 only emerges with increasing Res and is accompanied 
by an efficient and correlated momentum transport by the mean flow. Since the narrow 
maximum is of larger amplitude for Res = 2 x 10^, our simulations suggest that it will 
dominate at even higher Res- 
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1. Introduction 

The motion of a fluid between two independently rotating concentric cylinders, de¬ 
noted as Taylor-Couette flow (TCF), serves as a fundamental model system to study 
the influence of rotation on turbulence. The two main control parameters of the system, 
the rotation rates of both cylinders, can be combined into two parameters that capture 
the key physical processes. The differential rotation of the cylinders results in a shear 
gradient driving the flow, and the average of inner and outer cylinder rotation (denoted 
as mean system rotation) defines the magnitude of rotational influences which can either 
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Figure 1. Schematic representations of the two flows, (a) TC system seen from above in the 
reference frame rotating with ilrf- Only one third of the azimuthal circumference is shown, 
reflecting the actual computational domain for 77 = 0.5. (b) RPCF, as obtained from (a) for 
vanishing curvature, i.e. when 77 —>■ 1 for a constant gap width d. While the spanwise z-coordinate 
points out of the paper in both systems, the wall-normal coordinates r and y are antiparallel. 


enhance or attenuate the fluid motion depending on the rotation direction. This de¬ 
composition into differential and mean rotation is depicted in figure [TJa). Depending on 
the system rotation, the lam inar state becomes unstable due to a centrifugal instability 
I Rayleigh] 1917 : Tavloil 1 923[l or without a linear in s tability via a subcri tical transition 
I Wendt 1933 : Tavlod 1936t Burin fc Czarnockil[201^ IPeguchi et a^.ll2014ll similar to the 
transition in other wall-bounded shear flows. Beyond the laminar regime different com¬ 
binat i ons of shear and syste m rotation result in a multitude of distinct flow states ( ColesI 


196,51 : Andereck et al 


speeds (jRavelet et al. 


1986) that continue to sha pe the turbulence up to high cylinder 


2010; Huisman et al ]|20ii. 


The influence of the curvature of the cylinder walls can be captured by the ratio of 
the radii of the inner and outer cylinders, ry = rijro, as a third control parameter. If 
the radius ratio rj falls considerably below 1 , the curvature influence on the fluid motion 
becomes important. Conversely, the system curvature becomes small when 77 approaches 
1 while keeping the gap width — Vi constant, since the cylinder radii diverge in this 
limit. Then, the Taylor-Couette (TC) system approximates plane Couette flow subject 
to spanwise system rotation, c.f. figure (Hb) , commonly called rotating plane Couette 
flow (RPCF). We here study the limiting process to vanishing curvature for turbulent 
TCF and its dependence on the system rotation. For this purpose, we performed direct 
numerical simulations (DNS) of TCF for seven values of rj ranging between 0.5 and 
0.99. For each radius ratio, three turbulent shear rates and various mean rotation states 
were realised. We complemented these TC simulations with analogous DNS of RPCF to 
compare directly to the curvature-free limit. 

An important physical quantity in turbulent TCF is the torque needed to drive the 
cylinder s, which is connected to a radial transport of angular mome ntum by the fluid 
motion ()Marcusl l l984HEckhardt et a/.ll200nl : IPubrulle fc Hersantll2002li . In the turbulent 
regime, the torque can be characterised by its scaling wit h the shear strength modu- 


lated by a function that depends on the system rotation ( Dubrulle et al. 20^. This 


rotation-dependence of the torque featur es a maximum for counter-rotating cylinders as 


measured and calculated for rj = 0.72 (IPa oletti fc Lathronl 120111 : lyan Cils et al\ 12011 


Brauckmann fc Eckhardt 2013a; Ostilla et al. 201,3f l and for 77 = 0.5 ( Merbold et al. 


2 OI 3 I 1 . The toroue-maximizing rotation state has been linked to the occurrence of in¬ 


termittent turbulence near the outer cylinder that decreases the momentum transport 
efficiency ( yan Gils et al. 20121 : Brauckmann fc Eckhardt! l20136h and results from the 
stabilisation of the outer flow region for counter-rotating cylinders. At the same time, 
the flow near the inner cylinder remains centrifugally unstable. Since the curyature that 
causes this difference in stability between the inner and outer region yanishes for 77 —>■ 1 , 
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we also expect the rotation influence on the torque to change. Evidence for such de¬ 
viations in behaviour bet ween cases with n < 0.9 and larger n have also been noted in 
experiments with rj Ri 0.9 ( Ravelet et al ] |2010t[o^stilla-M6nico et a/.ll2014cJ) . However, the 
rotation-dependence in TC systems with a radius ratio significantly larger than 0.9 has 
not been explored, yet, and represents one of the main objectives of the present study. 
Moreover, we will study how turbulent characteristics, such as torques, mean profiles, 
and fluctuations, change during the transition from low-curvature TCF to RPCF and 
how TCF turns into RPCF. 

In the context of the transition to turbulence, the RPCF-limit of the TC system 
has already been studied in an approximation where certain terms in the equations of 
motion have been neglected ( Nagatal[l986L 119971 IPemav et al. 1992h . Nagata found that 
wavy vortices, which develop as an instability of streamwise invariant Taylor vortices, 
can be continued to non-rotating plane Couette flow where they form finite-amplitude 
solutio ns. Such wavy roll cel ls were also observed experimentally in RPCF at low shear 


rates (Hiwatashi et al. 


_ 2007f) . Moreover, by considering the full equations of motion, 

Faisst fc Eckhardt ( 2000ll were able to continuously study the transition from TC to 
plane Couette flow and found that already for radius r atios n below 1 some ch aracteristics 
of plane Couette flow can be observed. Furthermore. iDubrulle et ol ( 2005ll noted that 


also the linear stability criterion for the laminar state in RPCF (iLezius fc Johnston 
1976[l is approximated by the stability boundary of TCF ( Esser fc Grossmann 1996ll 
when the radius ratio rj tends to 1. Beyond the first bifurcati on towards turbulence , 
further transitio ns create a rich set of flow states also in RPCF (Tsukahara et al. 201(11 : 
Dalv et alll2014 1. reminiscent of the states described bv lAndereck et al. (1986) in TCF. In 
an extension of these bifurcation studies at low Reynolds numbers that revealed a number 
of similarities between the two systems, we here investigate the continuous transition of 
TCF to RPCF with a focus on turbulent characteristics at high shear rates. 

In the turbulent regime, the effect of rotation depends on the relative orientation of 
rotation and vorticity in plane Couette flow. Cyclonic rotation, i.e. when the rotation 
vector is parallel to the v orticity of the laminar ba se flow, was found to have a stabilising 
effect on the turbulence ( Komminaho et 0^ 1996 ) and to g enerate a striped pattern of 
coexisting laminar and turbulent regions (lTsukaharall2011ll . Similar turbulent-laminar 
pattern s were also observed at much higher shear rates by increasing the cyclonic ro¬ 
tation (jBrethouwer et aLN2012ll . On the contrary, anticyclonic spanwise rotation, where 
the rotation is antiparallel to the base-flow vorticity, has a de stabilizing; effect and was 
found to drive vortices reminiscent of turbu lent Taylor vortices (IBech fc Anderssonlll997 ; 
Barri fc Andersson 20I0t TsukaharrJbOIlll . For strong anticyclonic rotations, these vor¬ 
tices break up again and become disorganised. The influence of rotation on the strength 
of vortices is of substantial interest also in the present study since such vortices can 
enhance the momentum tra nsport and, therefore, cont r ibute t o the formation of torque 
maxima as demonstrated by iBrauckmann fc Eckhardt ( 2013611 . 

Investigations of RPCF for anticyclonic rotation furthermore revealed that the mean 
downstream velocity profile exhibits a linear section in the centre of the gap. In this 
region, which can extend over a considerable fraction of the gap, the v orticity of the mean 
profile compensates the impose d vorticity of the spanwise rotation (jBech fc Andersson 


1997; Barri fc Andersson 20inll . which means that the gradient of the central profile 


increases linearly with the imposed rotation rate. The arrangement of the profile to 
approximately realise zero absolute vorticity was recently also observed in a non-turbulent 
RPCF experiment, t hus demonstrating the generality of this behaviour for unstable flows 
( Survadi et al ] l20l4 . A physical explanation for these central profile shapes comes from 
the observation that mean profiles of zero absolute vorticity comply with neutral stability 
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in RPCF ( Barri fc Andersson 20101: Survadi et a,l. 2014 ). Such prominent changes of the 
profil e gradient with varying system rotation were also observed for the angular velocity in 
TCF ( Wendt 1933 : Ostilla et al. 2013 : Ostilla-Monico et al. 2014a). However, it remains 
unclear how the mean profiles in TCF compare to those in RPCF and whether they 
converge in the limit of vanishing curvature (ry —> 1). 

The purpose of the present study is to examine the continuous transition from TCF 
to RPCF and to compare turbulent characteristics between both systems taking into 
account the following aspects. First, we quantify the effect of the vanishing curvature, 
when rj tends to unity, by the decreasing strength of intermittent fluctuations near the 
outer cylinder and by the local stability for counter-rotating cylinders. Second, we show 
that as these intermittent fluctuations disappear, the rotation-dependence of the torque 
also changes. Third, we look for flow characteristics that prove the convergence of TCF 
to RPCF and seek the ry range in which the convergence can be observed. Fourth, the 
magnitude of the momentum transport by coherent vortices is analysed to assess their 
importance for the formation of torque maxima. In addition, we evaluate the effect of 
turbulent fluctuations on the momentum transport. Finally, we look for a convergence of 
streamwise velocity profiles and introduce a unified description of central profile gradients 
in TCF and RPCF. 

The paper is structured as follows. In ^we describe the unified framework in which 
we study TCF and RPCF, including the equations of motion, common control param¬ 
eters and their ranges, the momentum transport equations and numerical aspects. The 
vanishing of curvature effects is studied in followed in S]by an analysis of the total mo¬ 
mentum transport as well as the vortex-induced momentum transport. Finally, we discuss 
the rotation-dependence of velocity profiles and of turbulent fluctuations contributing to 
the momentum transport in ^ We close with a summary and some concluding remarks. 


2. Definitions and numerical methods 

Since there are several different combinations of parameters in TCF, we begin by 
discussing the limiting case of RPCF, so that we can then chose parameters that approach 
the ones in RPCF. As usual, RPCF refers to the flow between two parallel walls that 
move in opposite directions and are subject to system rotation in the spanwise direction. 
We describe this system in Cartesian coordinates {x, y, z) with the velocity field u = 
(ux, Uy,Uz) and the walls at y = —h and y = h moving in the x-direction with velocities 
—Uw and Uuj, respectively. The entire system rotates with angular velocity O = Hrfe^ 
around the spanwise axis as depicted in figure [Il[b). We chose d = 2h and Uq = 211^ as 
characteristic scales for all lengths and velocities. In these units and in a reference frame 
rotating with the Navier-Stokes equation and incompressibility condition become 

dtu + {u •V)u =-Vp - Raez X u +-^Au, V • m = 0, (2.1) 

Res 

where p denotes the non-dimensionalised pressure; the centrifugal force has been absorbed 
into the pressure. The motion is characterised by two dimensionless parameters, the shear 
Reynolds number Res and the rotation number Rq which measures the ratio between 
system rotation and shear. 


Res 


Uq d 

5 

V 


Rq. 


2 f2rf d 


( 2 . 2 ) 
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where v denotes the kinematic viscosity of the fluid. In contrast to the often used plane 
Couette flow Reynolds number Re = Uwh/v, the definition of the shear Reynolds num¬ 
bers Res is based on the full velocity difference and wall distance, so that Res = 4:Re. 

TCF is the motion of a fluid between independently rotating concentric cylinders. 
Its geometric parameters are the inner and outer cylinder radii and Tq and the axial 
height Lz- The radius ratio, defined as rj = r^jra, determines the curvature; it is thus the 
parameter that controls the transition to RPCF (which emerges for rj —>■ 1). The inner 
and outer cylinders rotate with the angu lar velocities LUj and uJn, respectively, and their 


ratio 11 = cVnltVj or its neg ative, a = —yt, ( van Gils et aLir2011 : Brauckmann fc Eckhardt 
2013Qt[Ostilla et ^1201311 is often used to define the rotation state of the s ystem. In order 
to stu dy TCF in the same framework as RPCF, we follow the analysis of lDubriille et al] 
(200^ and describe the flow in a reference frame that rotates with the angular velocity 
rirf. We make use of the cylindrical symmetry and introduce coordinates {r,ip,z) with 
velocities u = {ur,u^,Uz)- In analogy to RPCF, we choose Ilrf such that the cylinder 
walls move symmetrically in the opposite direction in the rotating reference frame, i.e. 

is fulfilled as depicted in figure [TKa). 


such that the condition u^{ri) = —Utp{ro) = 


This gives 


rfluii Ilrf) — ^o(^o f^rf) — 


rtUJt roUJo 


(2.3) 


Furthermore, we choose the velocity difference between the cylinder walls as the charac¬ 
teristic velocity scale Uo in the rotating reference frame, i.e. 


Uo = u^{ri) - u^(ro) = riCOi - roiOo + {ro - ri)r 2 rf 
2 ^ ^ 

= -r-- - WoUJo) ■ 

I rj 


(2.4) 


The last step in (12.41) results from substituting flj-f from equation (12.31) . With a view 
towards the variables in RPCF, we take the gap width d = r^ — rj as the cha racteristic 
scale to measure all lengths. Instead of the radius ratio n. lDubrulle et al\ (|2005ll specify a 
typical radius r = y/virp and define the curvature number Rc = d/r to characterise the 
curvature of the system. (The relation to rj is given in equation (j2.8l) below, where the 
definitions of all parame ters are collected.) Using these units, the equations of motion in 
cylindrical coordinates (|Landau fc LifshitallQSTT) transformed to the rotating reference 
frame can be written as 

dtu {u • V)m = — Vp — i?n 62 X M -I- —Am 


+i?c ■ 
Rc 




Res 


Uj‘UinG(n -(- 


1 


Res 


drU -I- -d^{ez X u) 


2 za 




Res 

~ r 

V • M = —Rc -Ur 
r 

with the modified Nabla and Laplace operators 

1 


V = Grdr 


-d^-\-ezdz, A = dr 




di 


(2.5) 


( 2 . 6 ) 


(2.7) 


Again the centrifugal force term is absorbed in a modified pressure p. The terms in the 
equations (12.51) -- (12.71) are arranged in a form adapted from Faisst fc Eckhardt ( 2000l) in 
order to clarify that in the limit ry —1 the equations of motion (12.5p and (12.61) indeed 
converge to the corresponding ones for RPCF (12.11) since the additional terms on the right 
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Figure 2. Transformation between the ratio of angular velocities = ujojuji and the rotation 
number Rn for various radius ratios rj. For TC systems with vanishing curvature ( 17 —^ 1), most 
71 -values result in almost no system rotation {Rn ~ 0 ) and Rn > 0 is achieved only for fj, close 
to one. The symbols mark the values of 71 and Rn realised by our simulations for Res = 2 x 10^. 


hand side vanish as Rq —>■ 0. The transition between Cartesian coordinates in RPCF and 
cylindrical coordinates in TC requires the identifications x = rip and y ^ r. 

The equations of motion show that th e TC system is characterised by three dimen¬ 
sionless parameters ( Dubrulle et al. 2005[l . 


Res = 


Uod 


2 ^ n \ r> 2nj.f d Rci -\- RCo 

{Rei-r]Reo), Rn =-p-r — = [^-V)- 


1 + 77 


Uo 


Rci — yRco 


Rc = - = — 

r 0? 


( 2 . 8 ) 


that result from (12.3(1 and (1^ . where Rci = TiOJid/v and Rco = roUJod/u denote the 
traditional Reynolds numbers that measure the dimensionless velocity of the inner and 
outer cylinders in the laboratory frame of reference. Note that Res and Rn are defined 
in strict analogy to RPCF (12.21) so that Rq, ha s the o p posite sign of t h e rota tion number 
used bv iDubrulle et al\ ()2005ll : iRavelet et al. ( 20ld) : IPaoletti et al\ ( 2012ll . Finally, in 
these units the laminar angular velocity profile in the laboratory frame given by the 
circular Couette flow reads 


Wlam(7’) — 2 


Rn ~ 1 ^ 


(2.9) 


with r ranging between and Tq. 


2.1. Parameters for the limit of vanishing curvature 
We are interested in the transition from TCF to RPCF, in the limit when the radius 
ratio rj = Vi/ro approaches one. This limit is often called the small-gap limit, since it can 
be obtained for fixed outer radius when the radius of the inner cylinder approaches that 
of the outer one and the gap width d = Vo — ri ultimately vanishes. However, in units of 
the width d, the radii increase like ri,ro 1/(1 — rj) and di verge as n approaches one : 


this is why this limit was called the case of “large radii” in Faisst fc Eckhardt ( 2n0(lll . 
Moreover, the curvature vanishes {Rc —>-0) as 77 —>■ 1, suggesting that curvature effects 
become less important. However, rotation effects can be maintained. 

For every radius ratio 77, the traditional Reynolds numbers can be expressed in terms 
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of the shear Reynolds number and the rotation number as 

^ R^o = -^ + Rn ( 2 . 10 ) 


and more specifically as Rci = Resl‘2 and Rco = —Res 12 in the non-rotating case 
(i?n = 0) for all r]. In contrast, for finite rotation {Rq ^ 0) the cylinder Reynolds numbers 
Rci and Rco diverge in the vanishing-curvature limit as 77 —>■ 1. As a consequence, the 
ratio of angular velocities, 


T = V 


RCq 

Rci 


- 77(1 - 77 ) -b ? 7 Rn 
(1 - 77 ) -I- pRn 


( 2 . 11 ) 


converges to 77 = —1 for i?n = 0 and to /r = 1 for any Rq 0 in the vanishing- 
curvature limit, as shown in figure [5] These limiting values illustrate that the transition 
between TCF and RPCF is singular in the traditional parameters. Ret, Rco, and /r, but 
not in the parameters Res and Rq that we will generally use here. We note that one 
exception to this parameter choice occurs in ^ where the effect of curvature for counter¬ 
rotating cylinders can be characterised more adequately by the ratio of angular velocities 
fi. Furthermore, while Rc measures the curvature and thus clarifies its disappearance 
(Rc —>■ 0) for 77 —>■ 1, we will identify the different TC geometries by the radius ratio 77 
as this parameter is more commonly used. 


2.2. Momentum transport in Taylor-Couette and plane Couette flow 


To study the mean turbulent characteristics of the flow, we calculate wall-normal pro¬ 
files of the velocity and other flow fields by averaging over surfaces parallel to the wall. 
The averages are denoted by {■■■). The surfaces are concentric cylinders in the case of 
TCF, (• • •) (r) = (• • and planes parallel to the bounding walls in the case of RPCF, 
(■ ■ ■) (y) = (■ ■ ')xz- Equations for the averages follow from the Navier-Stokes equation. 
The average of the 175 -component of the equation of motion (12.5p of TCF results in the 
continuity equation for the specific angular momentum C = rUip 

dt {C)r~^dr{rJc) = Q with Jc = (urC) — Rcg^r'^dr (ui), ( 2 . 12 ) 


where Jc denotes the angular-momentum flux in dimensionless units ( MarcusI Il984ll , 
u! = Upi/r the angular velocity and r~^dr{rJc) the divergence in radial direction. In the 
statistically stationary sta te, where dt (C) = O' th e transverse current = rJc becomes 
independent of the radius ( Eckhardt et al\\2QQ7lh . In physical units, the torque T needed 
to drive the cylinders, can be calculated from the dimensionless angular-momentum flux 
Jc sts T = pfiy^dRe^ A(r) Jc, with the dimensionless cylindrical surface area A{r) = 
2nrLz and the mass density of the fluid pf. The remaining prefactors compensate for the 
non-dimensionalization of (12.121) . We define the dimensionless torque per axial length. 


G = 


T 

2nLzPfi'^d 


RelrJc = Rel J^, 


and introduce the Nusselt number. 


(2.13) 


NU^ = G/Glan. = J^/Jflz,^, 


(2.14) 


that measures the torque in u nits of the torque of the laminar profile Giam l|Dubrulle fc Hersant 
2002t Eckhardt et a7ll2007(J) . 


Equations corresponding to (I2.12l) - (l2.14p can also be derived for RPCF. Averaging the 
downstream component of the equation of motion ( 12.111 results in a continuity equation 
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for the specific momentum in the x-direction 


dt{ux) - dyJ^ = 0 with = {uyUx) - Reg^dy{ux) 


(2.15) 

(see also Pooel 2005 §7-1), where J“ denotes the momentum flux in dimensionless units. 
It is independent of the wall-normal coordinate y in the statistically stationary state. 
In analogy to the torque T, the force F needed to keep the plane walls at a constant 
velocity can be calculated from the momentum flux J“ multiplied by the dimensionless 
surface area = L^Lz through which the flux passes, F = pfv^Re^g A^z Finally, 
we define the Nusselt number for the momentum transport in RPCF as in (I2.14L 


Nuu = F/Fi,^ = J“/J, 


lam 5 


(2.16) 


that measures the force in units of the value Flam for the laminar flow. 

This completes our set of relations between TCF and RPCF that we want to study here. 
We have identified correspondences between the continuity equations (12.121) and (12.151) . 
the torque and the force (T ^ F), the angular-momentum flux and the momentum 
flux {Jc O T“), the Nusselt numbers {Nuuj O Nuu), and the mean profiles that enter 
the diffusive part of the flux in (12.121) and (12.151) ((w) (r) o (ux) (y))- Since we study the 
statistically stationary cases only, we also average in time to improve the statistics. Using 
the dimensionless shear Reynolds number Reg and rotation number Rq (12.21) . we can 
now study the momentum transport in both systems within the same framework and in 
a non-singular manner. 

We conclude this section by pointing out a peculiar property of the relations that 
we will take up again in later sections. While the equations of motion m and (1^ 
depend on both driving parameters Reg and Rq, the momentum transport equations 
(12.121) and (12.151) do not explicitly depend on the rotation num ber Rq. Note that the 
same observation applies to RPCF ( Salewski fc Eckhardtl 1201511 . The rotation only in¬ 
directly influences the momentum fluxes Jc and J“ by changing the turbulent velocity 
correlations {(urC) and (uyUx)) and the mean profiles ((w) and (ux)), as we will discuss 
for the correlations in S (14.1.1I and 15.31 and for the profiles in 


2.3. Investigated parameter range 

Studies for different rj (or Rc) give information on when the properties of RPCF emerge 
from TCF, and studies of different Reg and i?n give information on the impact of shear 
and rotation on the turbulence. The range of parameters explored is indicated in Fig. [31 
They were selected as follows. 

For each curvature number Rc the evolution of the rotation dependence is anal¬ 
ysed with increasing shear, by realising various mean rotation states for three shear 
Reynolds numbers Res = 5 x 10^, 1 x 10^ and 2 x 10*^. Since we a dopt the results for 
■q = 0.71 from our previous study ( Brauckmann fc Eckhardt 2013ol ). the simulated val¬ 
ues of Res = 3899, 8772, 19737 deviate slightly from t he target values, as indicate d in 
figure |3Ka). Similarly, the RPCF results adopted from iSalewski fc Eckhardtl ( 2015 ) re¬ 
alise Res = 5.2 X 10^, 1.04 x 10^ and 2 x 10^. In the following, we will skip over these 
small differences in Reg, as they are not significant within statistical fluctuations. Our 
highest value of Reg lies in the range where the transition to the full y turbulent state was 
observed in experiment s for ry ^ 0.7 and stationary outer cylinder ([Lathrop et a/.l 11992 : 
[Lewis fc Swinn^ 19991) . We performed simulations for the range of rotation numbers 
Rq, where the occurrence of a torque maximum was observed for medium radius ratios n 
( Dubrulle et al. 2005t ^qletti fc Lathropj 2011 : van Gils et al. 201lL Paoletti et al. 2012 : 


Brauckmann fc Eckhardt] 2013ol : Merbold et al. 2013 ; Ostilla et al. 20131) : this allowed us 

to extend these studies to q close to one. Moreover, rotation numbers that correspond 
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Figure 3. Parameter space {Rc, Res, Rn) explored by DNS of TCF to study the limit of 
vanishing curvature, i.e. Rc —>■ 0 equivalent to g ^ 1. (a) For the shear Reynolds numbers 
Res and curvature numbers Rc indicated here, several simulations with Ra G (—0.3, 0.9) were 
performed, (b) The range of rotation numbers Rn explored is shown here for Res = 2 x 10“* and 
different curvature numbers Rc = (1 — g)/^/g- The curvature numbers correspond to the radius 
ratios g = 0.99, 0.98, 0.95, 0.9, 0.8, 0.71 and 0.5. The extra columns for Rc = 0 summarise our 
RPCF simulations. 


to strong counter-rotation fi < 0 were selected in order to study the fate of the tur¬ 
bulent intermittencY in the oute r region ( Coughlin fc MarcusI 199 (tI : van Gils et o,l. 2012 ; 
Brauckmann fc Eckhardtll2013Ml . Thus, the simulated rotation numbers lie mainly in the 
range —0.1 ^ ^ 0.6 as shown in figure[3Kb), except for Rc — 0.02 {g = 0.98) where we 

focus on a narrow range in i?n only for a study of the intermittency for counter-rotation. 
The case of no system rotation, i.e. Rq = 0, corresponds to counter-rotating cylinders 
with uJo = —goJi- 

For the highest shear Reynolds number Res = 2 x 10"*, the evolution of the rotation 
dependence is studied for seven curvature numbers Rc > 0 (TCF) and for Rc = 0 
(RPCF), as listed in table [TJ This investigation of the curvature effects compares tur¬ 
bulent characteristics, such as torques, mean profiles, and velocity fluctuations, between 
the different Rc studied. 

2.4. Numerical codes and tests of convergence 
For our investigations of TCF and RPCF we perform DNS that approximate a solu¬ 
tion to the incompressible Navier-Stokes equation and m, respectively. In 

these simulations we introduce a spanwise (axial) periodicity of length = 2 which 
is large enough to accommodate one pair of counter-rotating Taylor vortices. In addi¬ 
tion, a streamwise periodicity of length = 2n was assumed in the RPCF simulations 
while TCF is naturally periodic in this direction. However, we only simulated a domain 
of reduced azimuthal length that repeats ngym times to fill the entire cylinder cir¬ 
cumference. For g = 0.5 and g = 0.71 we tested that the reduced azimuthal length 
does not bias the computed torques for a stationary outer cylinder. The length Lj, was 
further increased for larger radiu s ratios, see table [H As dis cussed in more detail in 
Brauckmann fc Eckhardt (2013a); Ostilla-Monico et all ( 2015ll . the effect of these geo¬ 
metrical constraints of the domain are small for single point quantities, like the torque 
and profiles studied here. 
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V 

Rc 

n.sym 

Lx ^ L 

y 

X Lz 

Mx X My X Mz 

10'‘5, 

lO'^dy 

lO'* d, 

0.5 

0.707 

3 

7t 

X 

1 

X 

2 

62 X 

70 X 

190 

0.33 

0.74 

0.23 

0.71 

0.344 

6 

3.09 

X 

1 

X 

2 

78 X 

70 X 

110 

0.18 

0.54 

0.90 

0.80 

0.224 

7 

4.04 

X 

1 

X 

2 

94 X 

70 X 

94 

0.20 

0.33 

1.14 

0.90 

0.105 

15 

3.98 

X 

1 

X 

2 

94 X 

70 X 

94 

0.23 

0.37 

0.52 

0.95 

0.051 

30 

4.08 

X 

1 

X 

2 

94 X 

70 X 

94 

0.22 

0.35 

0.44 

0.98 

0.020 

78 

3.99 

X 

1 

X 

2 

94 X 

70 X 

94 

0.14 

0.20 

0.36 

0.99 

0.010 

99 

6.31 

X 

1 

X 

2 

158 X 

70 X 

94 

0.13 

0.30 

0.33 

1 

0 

- 

27t 

X 

1 

X 

2 

159 X 

76 X 

95 

0.42 

0.42 

0.38 


Table 1. Size of the simulations domain at the highest Res = 2 x 10'^ for the different radius 
ratios rj and corresponding curvature numbers Rc- L^, Ly and denote the streamwise (az¬ 
imuthal), wall-normal (radial) and spanwise (axial) length, and M^, My, give the highest 
order of spectral modes in the corresponding direction. The azimuthal periodicity Usym was 
introduced for the TC simulations to reduce the natural azimuthal length to = 27rra/nsyin 
with the mean radius Va = (rt + ro)l2. The amplitudes of the highest spectral modes dx, dy and 
Uz estimate the truncation error of the expansion in the corresponding direction. We listed the 
maximal amplitude values over all Rn. 


For the T C simulations w e emp loy a numerical code developed by Marc Avila and 
described in Meseguer et al. ( 20071 ). To accomplish large rotation numbers Rn ^ 0.6 
for radius ratios close to 1, both cylinders have to co-rotate rapidly and reach individual 
Reynolds numbers Rci and Rco of up to 6 x 10® while maintaining the shear Res = 2 x 10"^. 
Since these high absolute velocities introduce numerical instabilities, we perform these 
simulation cases in a reference frame rotating with t he angular velocity flrf defined i n 
(|2.3I) . For the RPCF simulations we use Channelflow ( Gibson et al. 2008t Gibson 2012 ). 
a pseudospectral code developed by John F. Gibson that was modified to include the 
Coriolis force, —Rne^ x u. Both numerical codes use an expansion of the velocity field in 
Fourier modes in the two periodic directions and in Chebyshev polynomials in the wall- 
normal direction, and employ a semi-implicit scheme for the time-stepping. We denote the 
highest order of modes employed in the streamwise, wall-normal, and spanwise direction 
by Mj;, My, and M^, respectively. 

The truncation error of the expansions can be assessed by the amplitudes of the highest 
modes dx, 5,, and Oz which are normalised by the globally strongest mode. Our resolution 
tests in Brauckmann fc Eckhardtl ( 2013al) revealed that these amplitudes have to drop to 
^ 10““^ in order to reach a converged torque computation. All present simulations fulfil 
this criterion. It turned out that the resolution in the axial direction could be reduced 
for larger ry while maintaining dz ^ 10““^. This seems plausible si nce the boundary layer 
for small t] is thinner at the inner cylinder than at the outer one (jEckhardt et oi.ll2n07Ml 
and thus requires a higher resolution. Nevertheless, we did not reduce the resolution in 
the radial direction for lar ge n. Moreover, our simula t ions sa tisfy two additional conver¬ 
gence criteria described by Brauckmann fc Eckhardt ( 2013ol l: Agreement of the torques 
computed at the inner and outer cylinder to within 0.5% and agreement of the energy 
dissipation estimated from the torque and from the volume averaged energy dissipation 
rate to within 1%. 


3. Curvature effects: Radial flow partitioning 

The first linear instability of TGF with counter-rotating cylinders is known to develop 
only in an annular region close to the inner cylinder while the outer region remains stable 
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( Tavloij[l923h . Inviscid calculations predict that the boundary between these two regions 
is given by the ne utral surface at the ra dius where the laminar Couette profile passes 
thr ough M m = 0 (IChandrasekhaiill96lL pp. 276-77). However, already the observations 
bv iTavloij (Il92,lh as well as viscous calculations sug gest that the structures fr om the 
inner partition protrude beyond r„ by a factor gfnl. lfeser fc Grossmaim ( 1996ll found 
it to be in the range between 1.4 and 1.6. Beyond the first instabilities, the initially 
stable outer partition beco mes susceptible to turbulent bursts at higher Reg ^ 2 x 10 ^ 
( Coughlin fc MarcusI 1 199611 which give rise to a radial partitioning of the turbulent flow 
that persists for much higher Res'- the coexistence of permanent turbulence in the inner 
partition and interm ittent turbulence in the outer one was observed for Res 10 ® 
201211 . 


(van Gils et ai 


In Braiickmann fc Eckhardtl ( 201811 1 we argued that the intermittency in the outer 


partition can only occur when the extended inner unstable region does not cover the 
entire gap. This resulted in the prediction 


Tpiv) = -V 


2 (a — iflrj + 0^ — 1 
( 2 o — 1)77 + 1 


with 0 ( 77 ) = (1 — 77 ) 


I (1 + 77)3 

2(1 + 377 ) 




(3.1) 


for the critical rotation ratio, pLc, so that flow partitioning occurs for /i < /Xc < 0 . 
The intermittent dynamics in the outer partition comes with a reduction in momentum 
transport and the formation of a maximum in torque as a function of fi. However, we 
already noted in our previous study that the bursting behaviour in the outer partition 
is a curvature effect that will not appear for 77 1 and, consequently, (13.11) will become 

invalid in this limit. 

In order to study what happens to the intermittency and the torque maxima, we 
characterize the stronger intermittent fluctuations in the outer partition by the standard 
deviation ac of the temporal torque fluctuations at the inner and outer cylinder for 
Res = 2 X 10^. This previously enabled us to iden tify the onset of enhanced oute r 
fluctuations as a function of p, for 77 = 0.5 and 0.71 ( Brauckmann fc Eckharddl20136ll . 
We here extend the analysis to 77 = 0.8 and 0.9. In the previous study, the relative 
fluctuation amplitude oqIG was assumed to be constant (i.e. independent of /i) at the 
inner cylinder. Therefore, the value of (JqIG at /r = 0 served as the reference value 
that the outer fluctuations exceed. However, since the data for 77 = 0.8 and 0.9 do not 
show constant inner fluctuations (figure|ll:,d), we here take the average inner fluctuations 
{uG IG) for strong counter-rotation as the base level. Figure [4] shows that the outer torque 
fluctuations exceed t his level when u decreases below t he critical value tJ-dv)- Similar to 
the procedure used in Brauckmann fc Eckhardtl ( 20136ll . we determine the critical values 


yidO-5) = -0.183 + 0.014, 
^idO-8) = -0.418 + 0.025, 


^ic{0■7l) = -0.344 + 0.050 
A7c(0.9) =-0.488 + 0.025 


(3.2) 


from the intersection of a linear ht to the outer fluctuations with the base level. For 
77 = 0.5 and 0.71, the re-evaluated onsets p-c conform with the old values within the un¬ 
certainties. Note that the difference between the fluctuation amplitudes is less pronounced 
for 77 = 0.9 than for 77 ^ 0.8 (cf. figure [4|), as expected from the weaker curvature in the 
first case. Moreover, the enhanced outer fluctuations due to turbulent bursting clearly 
differ from the situation with stationary outer cylinder or with co-rotating cylinders 
(/i ^ 0 ) where the fluctuation amplitudes are similar for both cylinders. 

Figure [5] shows that the prediction from (j3.ip agrees well with the onset of fluctuations 
in DNS data for Res = 2 x 10"^ (figure |4]). Moreover, Pp{r]) also agrees with the critical 
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Figure 4. Standard deviation aa (divided by the mean G) of temporal torque fluctuations calcu¬ 
lated at the outer cylinder (• O) and inner cylinder (A A) for a constant shear Res = 2.0 x 10"* 
and varying rotation ratio p. The dotted line marks the base level for the fluctuation at the 
inner cylinder for strong counter-rotation (average of open triangles). The solid line is a linear 
fit to 6 points in (a), to 5 points in (b), to 4 points in (c) and to 3 points in (d) which are 
marked as open circles. The vertical dashed line indicates the intersection point fic that defines 
the onset of enhanced outer fluctuations. 


rotation ratio for the occurrence of the bursting, dete c ted a s a bimodal distribution of 
angular velocities in the experiments of Ivan Gils et o,l\ (|2ni2l l for t] = 0.716. 

In figure UKa), we extend the analysis of the torque fluctuation amplitudes to the 
low-curvature radius ratios rj = 0.95, 0.98 and 0.99. For the sake of visual clarity when 
comparing all investigated rj we show the ratio between the standard deviation aa at 
the outer and at the inner cylinder. A value greater than one indicates enhanced outer 
fluctuations. Moreover, the rotation ratio ^ is shifted by the critical value in order 
to align the respective regions of enhanced outer fluctuations. The small difference for 
rj ^ 0.95 is the reason why we omitted a similar analysis as in figure |4] for the three 
highest radius ratios. To quantify the fluctuation asymmetry for strong counter-rotation, 
we dehne a typical fluctuation ratio ^ = ac{ro)/aG{ri) at the shifted rotation ratio 
^ — fj,c = —0.4 for each radius ratio. Figure [61(b) shows that this asymmetry measure. 
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V 

Figure 5. Onset of enhanced outer fluctuations. The critical rotation ratio fJ.c{v) for the onset 
is shown as a function of the radius ratio g. Numerical values for Res = 2.0 x 10^ are extracted 
from figure |4] and the value from the T®C-experiment (|van Gils et al.ll2012l l was deduced from 
a bimodal distribution of angular velocities. The line indicates the prediction fJ,p{g) from equa¬ 
tion m- As the radial partitioning disappears for 77 —>■ 1, we indicate the resulting uncertainty 
in the onset of enhanced fluctuations by continuing the curve with dashes for g > 0.9. 


(a) 


o 

b 


o 


O 

b 




Figure 6 . Enhanced outer fluctuations for Res = 2.0 x 10"^ as a function of the radius ratio 
g, assessed by the standard deviation of temporal torque fluctuations ac- (a) Ratio of a a 
computed at the outer and at the inner cylinder versus the shifted rotation ratio g, — fic, where 
Pc denotes the critical value for the onset of enhanced fluctuations from figure [5] For g ^ 0.95, 
Pp was subtracted instead of (b) Typical fluctuation ratio ^ = aciro)/adri) calculated by 
interpolation at g—pc = —0.4 marked by the vertical line in (a). The grey line represents a best fit 
to the data and indicates the disappearance of the fluctuation asymmetry near g^ = 0.99 ± 0.01. 


decreases monotonically with g towards a value of 1 which signifies equal fluctuation 
strengths near the inner and outer cylinder. The fluctuation asymmetry vanishes for a 
radius ratio of pc = 0.99 ± 0.01, which does not differ significantly from 1. This suggests 
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Figure 7. Radial dependence of the centrifugal instability for Res = 2 x 10"* and = —rj (i.e. 
Rn ~ 0) where the neutral surface is located in the centre of the cylinder gap. The real part of the 
large st eigenvalue A+ from the local stability analysis along Lagrangian path (lEckhardt fc Ya3 
Il995fl is shown for various radius ratios r/. The radial partitioning of stability decreases with 
increasing t; —^ 1 (indicated by the arrow). 


that small differences between inner and outer cylinder turbulence exist for all 77 < 1 . 
Therefore, in this case curvature has an effect for any 77 yf 1. 

The disappearance of the fluctuation asymmetry may also be related to the local 
stability properties of the flow. The results in figure [5] suggest that the partitioning into an 
unstable inner part and a stabilised outer part disappears as 77 —>■ 1. However, the neutral 
surface for counter-rotation, where the velocity profile passes through zero, discriminates 
the stability regions and exists for all 77 . Therefore, the radial partition ing has to vanish 
in a d ifferent way for 77 —>■ 1, as we will see from local stability results. lEckhardt fc Yaol 
( 1995h investigated the evolution of local perturbations to the laminar Couette profile 
wiam(u) along Langrangian trajectories. They calculated radially dependent eigenvalues 


A±(r) 


Res 


fci'(l + /3)± 


2wiam Q /" 2 \ 

- dr [r Wlamj ——5 

r ^ ' 1 + p 


1/2 


with 



(3.3) 


for the perturbation mo des with radial a nd axial wavenumbers fci and k^ (see also the 
local stability results by Dubrullel ( 19931) ). Note that the wiam-dependent factor in the 
second term of (I3.3|l corresponds to Rayleigh’s stability discriminant. In the following, we 
focus on the largest (most unstable) eigenvalue that can be calculated for given Res and 
p by maximizing A+(r) at the most unstable radial position r = Vi. For that purpose, 
we minimize the viscous damping (first term in (13.3p ) by selecting the smallest possible 
radial wavenumber 


ki 


nj{ro-ri), 
7T/(r„-ri), p<0 


with 



(3.4) 


Subsequently, we determine the wavenumber ratio /3max that maximizes equation ()3.3I1 . 
With these^optimal values for ki and /3, we study the radial variation of the largest 
eigenvalue A+. 

To analyse the local stability for strong counter-rotation, figure [7] shows radial profiles 
of A+ for Res = 2 x 10"* and the specific rotation ratio p = —rj. This value is comparable 
to the value p — pc = —0.4 used in the^study of the fluctuation asymmetry in figure [H 
The positive real part of the eigenvalue A+ inside the neutral surface indicates instability 
in the inner partition while the outer partition is vicously damped. Most importantly. 
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the radial variation of the local stability decreases with p which is reminiscent of the 
disappearance of the radial differences in the fluctuation behaviour. For large rj, the 
eigenvalues move in opposite directions: in the inner partition, they become less unstable 
and approach neutral stability from above, whereas in the outer partition, they approach 
the neutral value from below. 

The results of this section show that the radial partitioning of the flow that gives rise 
to enhanced outer fluctuations conforms with the prediction for 77 ^ 0.9. However, 
this curvature effect as well as the radial variation of the local stability disappear with 
vanishing curvature for 77 1. The disappearance of the intermittent bursts in the 

outer partition is also relevant for an understanding of the variation of torque, since 
the appearance of these bursts has been linked to th e emergence of a torque maximum 
( van Gils et aLl[201^ iBrauckmann fc Eckhardtll20136ll . The absence of this intermittency 
should therefore change the variation of torque with rotation, as we now discuss. 


4. Variation of the momentum transport 

We now turn to the rotation dependence of the mean torque with increasing Res and 
with variations in radius ratio. Figure | 8 ] shows the Nusselt number Nu^j, the torque 
in units of its laminar value. For the wide-gap TC system with 77 = 0.5, the torque 
maximum initially occurs at i?n 0.3 for Res = 5 x 10^ and t hen shifts to Ro =0.26 for 
Res = 2 X 10^ correspondin g to p, = —0.20 as analysed before ( Brauckmann fc Eckhardtl 


20136 : Merbold et a/.ll2013[) . For rj = 0.71 this shift in the maximum is more pronounced, 
with i?n varying from 0.27 at Res = 5 x 10 ^ to 0.11 at Res = 2 x lO'^ (figure [8b ). For 


even higher Res, the experiments of Paoletti fc Lathr^ (20111 and van Gils et aZI (|2nil 


1201 2II show no further shift in the position of the torque maximum so that Rq = 0.11 is 
likely to be close to the asymptotic value for high Res- As discussed in ^ these high-i?es 
torque maxima for 77 = 0.5 and 0.71 coincide with the onset of intermittent bursts in the 
outer partition. 

The situation changes drastically for low-curvature TCF with 77 ^ 0.9 as shown in 
figures [ 8 l[d)-(f): instead of a single maximum, one notes a broad torque maximum near 
= 0.2 and a second narrow maximum near i?n = 0.02. This narrow maximum 
emerges with increasing shear and, at Res = 2 x 10 ^, is similar in magnitude to the broa d 
maximum. Indications for the narrow maximum were first seen bv iRavelet et a/] ( 2010f tl 
in a low-curvature TC experiment (77 = 0.917): Their data show a very slight bump in the 
torque near Rq = 0.02 for Res = 1.4x 10 "* and 1.7x 10 “^ and a monotonic increase with Rq 
for Rq > 0.05. However, their figure 8 only shows torque measurements for Rq < 0.125 
so that the broad maximum found here at Rq = 0.2 lies outside their investigated range. 
The presence of two maxima, a br oad one and a narrow one th at emerges with increasing 
Res, was also observed in RPCF ( Salewski fc Eckhardt 2015ll . They compare well to the 
Nuui maxima in low-curvature TCF as demonstrated for 77 = 0.99 in figure [SDf). 

Finally, the torque for 77 = 0.8 (figure H}:) shows a behaviour intermediate between 
that of the narrow-gap and the wide-gap TC systems. While the broad torque maximum 
known from systems with larger radius ratios is present, the second maximum around 
Rq = 0.02 is not as narrow and as clearly visible as for 77 ^ 0.9. On the other hand, the 
extrapolation of equation which captures the maximum for 77 = 0.5 and 0.71, to 

the case 77 = 0.8 predicts a maximum at fj.p = —0.419 and Rn.p = 0.069. This is close to 
the low-curvature narrow maximum and suggests that the two will interfere. Indeed, the 


f The rotation number Ro employed by these authors is of opposite sign, i.e. Ro = —Rn- 
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Figure 8. Variation of the momentum transport {Nuui) with changing mean rotation (Rn) 
calculated in DNS for the radius ratios rj specified in the sub-plots. The simulations are performed 
for Res = 5 x 10® (■), Res = 10^ (•) and Res = 2 x lO'^ (A). The open symbols in (f) show 
results from simulations of RPCF. Note the increased i?n-scale in sub-figures (a) and (b). 
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Figure 9. Comparison of the momentum transport {Nuui) for various radius ratios as a function 
of the mean rotation (Rn) at Res = 2 x lO'^. The values are compen sated with a scaling in Res 
employing the scaling exponent measured by Ivan Gils et gp (120121 1. Sub-figure (b) details the 
narrow maximum in (a). For purposes of clarity, the results for g — 0.98 are only included in (b) 
where the ones for g — 0.5 and g = 0.71 are left out. The open stars labelled by ?7 = 1.00 show 
results from RPCF simulations. The inset in (a) com pares our DNS results for g = 0.9 to the 
T^C-experiment with g = 0.909 and Res of a few 10® llOstilla-Monico et aZ.ll^llofl and reveals 
a slight shift of the maximum with increasing Res- 


magnification in figure [Hfb) of the region in around the low-curvature narrow peak 
shows both, a slight maximum at Ra = 0.02 and a broader one around Rq, = 0.05. 

In figure [9l we show the torque as a function of the rotation number for various radius 
ratios. In order to co mpensate for the slight ly different Res of the compared torques, we 
use the observation bv iDubrulle et al\ ( 2005ll that the torque variation can be decomposed 
into a scaling with Res modulated by an amplitude that describes the dependence on 
Ra and we thus normalise the torques by For g ^ 0.9, the Nusselt numbers 

for different radius ratios collapse. The figure also contains data for RPCF which also 
agree nicel y with those of l ow-cu rvature TCF. Data collapse is also observed with other 
quantities. Dubrulle et all 1 2005t) previously noted that the torque G normalized by its 
value for stationary outer cylinder G{u = 0) varies only little with g when plotted as 
function of Rq. Similarly. IPaoletti et all ( 2012ll found a collapse of another torque ratio 
G/G{Rq = 0) as a function of Rq for various radius ratios. Our figure [5] shows that 
up to a scaling with Res a normalisation by the laminar torque suffices to compensate 
for the curvature dependence and to achieve the collapse for g ^ 0.9. Note that p, = 0 
corresponds to Rq = 1 — ry so that the normalisation by G{fj, = 0) has an additional 
77 -dependence compared to the normalisation by G{Rq = 0). The collapse of the data 
in all three cases is of similar quality and does not favour one particular normalisation. 
However, plots with respect to /i (not shown here) do not give such a collapse, confirming 
that Rq is the appropriate parameter in which to describe the transition from TCF to 
RPCF. 

For lower 77 , the collapse of the torques is limited to the region Rq > 0.1 for 77 = 0.8 
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Figure 10. Ratio of angular velocities of maximal torque pmax as a function of the curvature 
number Rc = (1 — 'ri)/^/rj. The open circles, triangles and squares summarise our numerical 
results for the detachment, broad and narrow maximum, respectively, at Res = 2 x 10"^ from 
figure [9l The maximum locations /imax were calculated from quadratic fits to the peaks in 
NuuiijJ,), except for the narrow maximum for r] ^ 0.95, which is not well approximated by a 
parabola and therefore determined by the location of the largest data point. The small circles 
indicate the experimental results bv lOstilla-Monico et al\ ( 2014af) for Res ~ 10^ to 10^. These 
are complemented by the predicti on Uv(ri) from equatio n (|3 l|l ijBrauckmann fc Eckhardtll20136f ) 
and by the angle bisector fib{v) bv ivan Gils et aOll2012ti . In addition, two for constant 

rotation number Rn = 0.2 and Rn ~ 0.02 are shown, corresponding to the broad and narrow 
torque maxima. 


and to Rq > 0.25 for 77 = 0.5 and 0.71. For these rotation numbers, turbulent bursting 
does not occur in the outer partition. For smaller Rq, where the bursting in the outer 
partition occurs, the Nusselt numbers in figure [Oja) depend on rj, suggesting that only 
this radial flow partitioning introduces the strong curvature dependence of the torques. 
In conclusion, the impact of the global rotation (parametrised by Rq) on the turbulent 
momentum transport becomes independent of rj if the flow is not partly stabilised by 
counter-rotating the outer cylinder or if this curvature effect becomes negligible for large 
radius ratios. 

The position of the torque maxima for Res = 2 x 10"^ and 77 ^ 0.9, which in Rq 
is given by constant rotation numbers Rq = 0.2 and Rq = 0.02, translates into 77 - 
dependent rotation ratios /x, which are represented by the thick solid lines in figure 
[TOl The torque maximum locations fimax from the DNS (triangles and squares) closely 
follow these lines for 77 ^ 0.9. In particular, this reveals that for 77 = 0.99 both the broad 
and the narrow torque maximum occur for co-rotating cylinders, i.e. fj, > 0. Therefore, 
the narrow maximum cannot be related to the detachment mechanism that causes the 
torque maximum at counter-rotation for 77 = 0.5 and 0.71. For the latter radius ratios, 
the maximum location /Xmax agrees with the predictive line fJ-pir]) in figure ITOl Finally, the 
intermediate radius ratio rj = 0.8 shows indications of all three torque maxima close to 
the lines in figure [101 with the detachment and the narrow maximum lying close together 
and being not as pronounced as the broad maximum, cf. figure |9l 

In addition, figure [101 shows the torque maxima identified for various radiu s ratios 
in experiments at much higher Reynolds numbers (jOstilla-Monico et 0,1.112014 oi l. While 
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the maximum positions /imax for p — 0.714, 0.769 and 0.833 conform with the trend of 
the detachment maximum, the one for rj = 0.909 clearly deviates from both /ip and our 
simulations; this will be discussed in detail at the end of th i s sect ion. Moreover, the angle 
bisector line /ib(r?), that was suggested by van Gils et all (2012) for the position of the 
maximum /tmax, lies close to the experimental t orque maxima. Howe ver, p,b{vi) clearly 
deviates from the torque maximum for rj — 0.5 (iMerbold et a/.l 120131) and differs in its 
functional behaviour for ?/ —>■ 1 from our DNS results. 

For the range of shear Reynolds numbers Res ^ 10"^ investigated here, the rotation 
dependence of the torque significantly changes with increasing Res, c.f. figure[5] It is likely 
that this transformation process continues beyond Res = 2 x 10"^. However, concerning 
the detachment torque maximum for 77 = 0.5 and 0.71 it was shown that, after an 
initial transformation with increasing shear, the rotation dependence of the simulated 
torques at Res = 2 x 10 ^ already compares well to experiniental torque measurenients a t 
much higher Res ~ 10^ to 10® ( Brauckmann fc Eckhardt 2013Qt Merbold et a/.|[2013h . 
On t he other hand, conc e rning t he two torque maxima in low-curvature TCF (ry > 
0.9b lOstilla-Monico et al. (2014a) find in their torque measurements for rj = 0.833 and 
0.909 that the power-law exponent 7 of the torque scaling Nu,^ ^ Res'^ depends on 
the system rotation for Res below a few 10 ®, which indicates that the transformation of 
the rotation dependence still continues up to Res ^ 10 ®. As regards the second narrow 
torque maximum, we expect that it grows further in relation to the broad maximum for 
i?es > 2 X 10^ and that this may shift the torque-maximisi ng rotation number a bit. 

The torque maximum from the experiment at ?/ = 0.909 (jOstilla-Monico et ai!.ll2014af ) 
supports this assumption as shown in the inset of figure |9l[a) which compares experiment 
and DNS: The experimental torque-maximising rotation ratio /tmax = —0.34 corresponds 
to i?a = 0.04 which only slightly deviates from Rq = 0.02 for the narrow maximum at 
lower Res = 2 x 10"^. In conclusion, apart from small changes with Res, our simulations 
already capture the beginning of the general turbulent behaviour at much higher shear 
rates. 


4.1. Torque due to vortical motion 

In this and the following sections, we investigate the momentum transport characteristics 
that underlie the rotation dependence of the Nusselt number discussed above. Vortical 
motions, such as Taylor vortices and their turbulent remnants, are known to effectively 
contribute to this transport by moving fast fluid from the inner cylinder outwards and 

slow outer fluid inwards. Their effects on the turbulence have been discussed previously _ 

( Lathrop et a/.l[l9^ ; Lewis fc Swinne\^ll999t Martfnez-Arias et a/.ll20l3) . iBrauckmann k, Eckhardt 
(|20136i) quantified the link of the vortical motion to the torque and found that the torque 
contribution of the mea n vortical motion was larges t near the torque maximum. Flow 
visualizations of DNS by Salewski fc Eckhardt ( 20151) show distinct vortical states which 
underlie th e narrow and broad m axima for RPCF. Furth ermore, using experiments a t 
Res ^ 10® ( Huisman et al\ 2014 ) and DNS at Res ^ 10® ( Ostilla-Monico et al. 20146 ). 
these structures have been detected over a range of rotation numbers where their presence 
is associated with a single torque maximum. 

To investigate the effects of vortical structures to the torque, we first extract the mean 
vortical motion underlying the turbulence in the full simulation and measure its contri¬ 
bution to the momentum transport. The mean vortices consist of temporally and stream- 
wise averaged turbulent Taylor vortices and are a nalogous to the “lar ge-scale circulation” 
found in turbulent Rayleigh-Benard convection ( Ahlers et aLl 120(391) . Next, since Taylor 
vortices do not depend on the streamwise direction, we also consider the torque which 
results from streamwise independent flow. We then summarize our results in figure [TTJ 
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4.1.1. Mean vortical motion 

We follow the procedure introduced in Brauckmann fc Eckhardtl ( 20136ll and decom¬ 


pose the flow of the full DNS into a mean contribution u = ^ (or u = (u)^ ^ for 


/ X,t 


RPCF) that includes the mean vortical motion in the rz-plane (yz-plane) and into the 
turbulent fluctuations around these mean vortices u" = u — u. Substituting this velocity 
decomposition into the momentum flux equation (I2.12p results in a separation of the 
torque into a contribution that is caused by the mean vortical motion G and a second 
contribution that is due to turbulent fluctuations G", \.e. G = G -\- G” with 


G = 


Rei 


C2 - n 


[ - Res^r^dr rdr, 

J VT ^ ^ 


G" = 


Rei 


/ n 
C ‘^2 


r2 - ri 






■ dr. 


(4.1) 


Analogous expressions for the decomposition of the driving force into F = F + F" can be 
obtained for RPCF by substituting the velocity decomposition into the momentum flux 
equation (12.151) . To capture the correct amplitude of the turbulent Taylor vortices, we 
hx their spanwise position during the temporal average of u and in (14.11) by correcting 
a potential spanwise drift. Such a drift would lead to a cancellation of the large-scale 
motion over time, even though it might be strong. Note that while the total transport G 
is constant over r, the individual terms in (14.111 vary with the radius. This motivates the 
additional radial average in dni in order to quantify the typical strength of the torque 
contributions. Since we are interested in the mean vortical motion which dominates in 
the central region and not in the boundary layers, we restrict this average to the range 
between (ri — ri)/d =1/4 and (r 2 — ri)ld = 3/4. 

In hgureim the torque due to the large-scale circulation G is compared to the total 
torque G. For all cases, G exhibits a maximum. For r] ^ 0.8 and RPCF, the maximum 
of G, where the momentum transport by the mean vortices is most effective, nearly 
coincides with the narrow maximum of the total torque G. Moreover, it is apparent that 
a distinct narrow peak also occurs in mean-vortex torque G for large 77 . This suggests 
that the narrow, high-i?es maximum is linked to an efficient large-scale circulation. For 
lower radius ratios, i.e. 77 = 0.5 and 0.71, the coincide nce between the maxima in G and 
G is not as significant; nevertheless, as the analysis in iBrauckmann fc Eckhardtl (j20136[ l 
demonstrates, this coincidence does exist for 77 = 0.5 and 0.71 when the radial average 
in (HID covers the complete radial gap instead of the central region. Finally, the mean- 
vortex torques G in figure [Tl] reveal that the mean Taylor vort ices occur and grow in 
strength for i?n > 0 as previously observed in the experiment by Ravelet et al ] (I 2 OIOII . 


4.1.2. Importance of streamwise invariant structures 

We perform simulations that force the flow to be invariant in the downstream direc¬ 
tion by taking no Fourier mode in the (/?—direction (x—direction for RPCF). For TCF 
this corresponds to axisymmetric simulations. Note that this procedure differs from 2D 
simulations in a cross-section because three velocity components are still active, however, 
allowing only for spatial variations in two directions (wall-normal and spanwise). These 
simulations of streamwise invariant vortical flow result in torques, denoted as G 2 D, that 
we compare to the total torques G in hgure [TTJ The torque G 2 D also shows a broad 
maximum, which has an 77 -dependent shift in its location: for 77 ^ 1 , the maximum is at 
Rq, = 0.18; and for 77 = 0.5, it is near Rq = 0.35. For some radius ratios, namely 77 = 0.8, 
77 = 0.9, this places the G 2 D-iiiaximum at nearly the same as the broad maximum 
in the total torque; however, the streamwise invariant simulations overestimate the am- 
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Figure 11. Comparison of the total torque G (■) to the torque that results from the streamwise 
invariant flow 6*20 (®) and to the torque due to the mean flow G (A) for Res = 2.0 x 10'*. All 
torques are measured in units of the laminar torque Giam giving Nusselt numbers Nui^. The 
open symbols in (f) show corresponding results from RPCF simulations. 
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plitude of the maximum. In fact, for rj ^ 0.8 and RPCF, the G 2 D-niaximum resembles 
the broad maximum in the total torque. Furthermore, it also appears to agree with the 
plateau that occurs in the total torque beside the actual torque maximum for ry = 0.71 
(figure fTTbl. This suggests that the plateau for Res = 2 x 10^ as well as the broad 
maximum at lower Res (cf. figure [Dd) follow from the same mechanism that causes the 
broad maximum for rj ^ 0.8. 

The analysis reveals several features of the contribution of vortical motion to the torque 
for different degrees of curvature. The streamwise-averaged mean flow (G) appears to 
reproduce the narrow torque maximum, but this is not reflected by the streamwise- 
invariant flow (G 2 d)- This seems to be counter-intuitive when considering that u is 
streamwise invariant like the simulations underlying G 2 D, but neither the detachment 
maximum in the total torque for rj ^ 0.8 nor the narrow maximum for rj ^ 0.8 occur in 
the torques from the streamwise invariant simulations. As a consequence, these maxima 
must arise from more complicated flows which allow for streamwise fluctuations and 
cause an additional increase of the momentum transport. The reasons for the strong 
momentum transport differ between the broad and narrow torque maximum as section 
>15.31 will illustrate: strong streamwise vortices as well as a high correlation between the 
radial flow u^. and the angular momentum C in these vortices are required for an effective 
momentum transport and, thus, high torque. 

5. Flow characterisation 

In section |T] we discussed that the momentum transport converges to a universal be¬ 
haviour for the low-curvature TCF {rj ^ 0.9), in that the rotation dependence becomes 
independent of rj. Furthermore, the dependence on Rq agrees with that observed in 
RPCF. In the following, we analyse how the convergence for rj ^ 0.9 extends to other 
flow characteristics such as mean profiles and turbulent fluctuations. Moreover, we com¬ 
pare our TC results to the behaviour in RPCF. 

5.1. Angular momentum profiles 

We analyse here the mean profiles of the specific angular momentum C = ru^p which is 
radially transported between the cylinders according to the continuity equation (12.121) . 
To compensate for the varying global system rotation, we show rescaled profiles C = 
((£) — Co) ! {Ci — Co) in figure [T^ where Ci and Co denote the specific angular momen¬ 
tum of the inner and outer cylinder. As a result, the profiles for various radius ratios rj 
that belong to the same rotation number Rq collapse as long as the flow is not subject 
to the radial partitioning of stability for moderate counter-rotation (cf. section [3]): The 
rotation number = 0.5 corresponds to co-rotation {ji ^ 0) for all investigated rj., and 
all profiles of C in figure IT^dl collapse well. However at Rn — 0.2, only the simulation for 
?7 = 0.5 with ji = —0.22 shows the bursting in the outer partition and a deviating profile 
in figure [r^ cb Similarly at i?n = 0.02 (figure [T2 )d), the profiles for ry < 0.9 corresponding 
to /r-values in the bursting regime clearly deviate from the profiles for rj = 0.95 and 0.99 
with ji > jip which are closer together. Finally, no collapse is observed for Rq = 0 which 
always corresponds to exact counter-rotation. In addition, figure [12] shows Ux profiles 
from RPCF with the system rotation subtracted, i.e. (ux) — Rnjj, that agree well with 
the collapsed C profiles from TCF when rescaled to the same interval. The reason for this 
correspondence between angular momentum profiles in TCF and {ux) — Rny in RPCF 
will be discussed in the next section. 

The C profiles tend towards a universal shape for increasing i?o, as the radial parti¬ 
tioning of stability disappears. The specific angular momentum in the centre becomes 
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Figure 12. Profiles of the rescaled specific angular momentum £ = {{£) — Co) / {Ci — Co) for 
the mean system rotations defined by Rn and for six different radius ratios. The corresponding 
RPCF profile {ux) — Rny rescaled to the interval [0,1] is shown by the dots. Each figure contains 
profiles for g = 0.5, 0.71, 0.8, 0.9, 0.95, and 0.99 as listed in (d). Some of the profiles collapse, 
especially in (c) and (d). In (a) and (b) the profile centre level decreases monotonically with 
increasing t; —^ 1. The dotted line shows the arithmetic mean angular momentum C — 0.5. 


nearly flat (with a slightly positive slope) and approaches C = 0.5 which corresponds 
to (C) = (Cj + Cn) !^ and indicates that the angular momentum is well mixed. Al- 


ready 


IWendtf 


(|1933f) noted that in the unstable regime the angular momentum be¬ 


comes almost constant in th e centre which was also repeat e dly observed for the case 
of stationary outer cyl i nder (IWattendorl 1935 : Tavlor 1935 : ISmith fc Townsend Il982l : 
Lewis fc Swinnev 19991 : Donellio^ 1. Flat £ profiles resemble the situation in Rayleigh- 


Benard c onvection where the transport of heat results in a well mixed temperat ure in 
the bulk ( Tilgner et al. 1993 ; Kerr 199fil : Brown fc AhlersI 2007 ; Ahlers et flZ.ll20l3 l . 


5.2. Mixing of angular momentum and neutral stability 

The observation that the specific angular momentum £ = ru^p becomes almost flat in 
the middle is remarkable and calls for an explanation. Moreover, it is of interest whether 
an analogous phenomenon can be found in RPCF. We will approach both questions 
by the stability properties of the mean profile around which the turbulence fluctuates. 
As long as the mean profile ((w) (r) in TCF and (ux) (y) in RPCF) is unstable due 
to a strong a gradient in the centre, it further enforces turbulence which increases the 
mixing of angular (or linear) momentum; the mixing flattens and therefore re-stabilises 
the profile. However, a stable, low-gradient profile decreases turbulent mixing, so that 
the profile becomes steeper and eventually unstable again. Therefore, we assume that 
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both processes balance in the statistically averaged sense and the mean turbulent profile 

becomes neutrally stable in the centre. __ 

For TCF with stationary outer cylinder. IWattend^ ( 1935 1 and Tavlor ( 1935ll already 
noted that the mean angular momentum profile in the middle tends to wards one of neu¬ 
tral stability according to Rayleigh’s inviscid criterion (lRavleighlll917ll , which predicts a 
constant angular momentum profile for this case. Neglecting viscous effects in the central 
part of the flow becomes justifiable at high Reynolds numbers. The flat C profiles in figure 
[12] indeed realise Rayleigh’s neutral stability criterion in the centre, thus, supporting the 
aforementioned assum ptions. More precisely, Rayleigh’s discriminant for neutral stability 
( Chandrasekhai][l96lh reads 


0= — dr = 


i?n — 1 -l- 


= i?n - 1 + 


-drC 


Rn-l 


rdr 


(5.1) 

(5.2) 


here expressed so that it serves as a condition for a neutrally stable turbulent mean profile 
in either £ or a; in the middle. For that purpose, we substituted the angular velocities 
representing the mean profile by the laminar solution uJiamir) from equation (1^ because 
in the centre wiam provides a reasonable approxi mation to the magnitude of the tu rbulent 
mean profile as evidenced by the comparison in Brauckmann &: EckhardO ( 2013a ). Since 
in the middle r is close to the geometr i c radi us r = y/riVo, equation (15.21) implies that 
rdrco ~ —Rn for Rq ^ 0. lOstilla et a,l. (2013) obtained a similar linear relationship by 
assuming that the Coriolis force term balances the convective term in the (^-component 
of the e quations of motion dSSj). Equations dUD and (15.211 then show that both force 


balance ( Ostilla et al. 2013h and inviscid stability ( Ravleigd 1917 ) predict a profile with 
flat angular momentum. 

The equivalent stability formulation in RPCF was derived from an a nalogy between 
rotating shear flows and buoyant flows ( Bradshaw 19691 : Tritton 1992ll : The stability 
expressed by a quantity similar to the Richardson number in stratified flows results in 
the neutral stability condition 


Rn (Rn - dyU^) = 0 


(5.3) 


which is approximately contained in (15.21) since r/r « 1 in the middle, and rdrOJ corre¬ 
sponds to —dyUx - The difference in sign originates from the fact that the radial coordinate 
in TCF is antiparallel to the wall-normal coordinate in RPCF, cf. figured] Also for RPCF 
we assume that the turbulent mean p rofile in the m i ddle l aecomes neutrally stable and 
fulfils (15.31) . as proposed previously by Survadi et al\ 1 2014h . In conclusion, the negative 
profile slope —dy (ux) in RPCF corresponds to the angular velocity gradient rdr (w) in 
TCF. Moreover, comparing (|5.3I) to EU) reveals that the combination (i?n — dy (ux)) 
in RPCF plays the role of the angular momentum gradient r~^dr (£) in TCF, and we 
therefore compare the modified RPCF profile {ux) — Rny to the angular momentum 
profiles in figure [12] 

To test the assumption of neutrally stable mean profiles, the dependence of the profile 
gradients on Rq, and the proposed analogies between TCF and RPCF, we evaluate the 
relevant profile gradients in the central region for Res = 2 x 10^ and various radius 
ratios. Figure fTSl a) shows that the specific angular momentum becomes almost fiat with 
a slightly positive slope for Rq > 0.02, except for rj = 0.5 which is well mixed only for 
Rn ^0.1. Experimen ts show that the slightly positive angular moment um slope persists 
for much higher Res ( Smith fc Townsend 1982 : Lewis fc Swinnev 1999[) . In addition, for 
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Figure 13. Profile gradients versus the rotation number Rn for Res = 2.0 x 10^ and various 
radius ratios g. (a) radial gradient of the specific angular momentum profile {£,) in TCF and 
the corresponding quantity {Rn — dy (ux)) for RPCF, (b) gradient of the angular velocity profile 
(oj) in TCF and the negative downstream velocity gradient —dy (ux) in RPCF, (c) detail of the 
maximum in (b). The black dashed line in (b) indicates the estimate rdr (lo) = —dy (ux) = —Rn 
from the neutral stability conditions (1^ and (lOll . All gradients are averaged in the central 
region (r — ri)/d G [0.4, 0.6]. 


most i?n > 0, the profile slopes measured by r~^dr (£) collapse for different rj, which 
demonstrates the universality of the weak i?o-dependence. Furthermore, the quantity 
{Rn — dy (ux)) from RPCF also shows this high level of agreement expected from the 
neutral stability analogy (EH) and (15.31) . In contrast, the profile gradients for Rq < 0 
{Rn < 0.1 for ?; = 0.5) clearly vary with the radius ratio. For these rotation numbers, 
a pronounced stable partition exists which introduces a strong curvature dependence as 
discussed for the torques in section g) 

The flattening of the specific angular momentum implies that the a ngular velocity gra¬ 
dients d epend on Rn as shown in fjgure fT^T b.cl and in similar figures in lOstilla-Monico et al. 
( 2014cJl . In the unstable regime with Rn > 0, the profile slopes rdr (w) approach the 
neutral stability estimate rdr (w) = —dy {ux) = —Rn as ry tends to 1. Moreover for low- 
curvature TCF {rj > 0.9), they coincide with the Ua,-profile gradients in RPCF which sup¬ 
ports the proposed analogy between (w) a nd jufl)- The approximate s calin g {ufl) ^ Rn 
in RPCF conforms with observations bv iBech fc AnderssonI (1993) and Survadi et a,I 
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( 20141 ). However for smaller 77, the (a;)-gradients depend more strongly on the curvature 
than the (£)-gradients, consistent with the additional radial dependence in condition 
(ESI). Furthermore, we note that a slight local super-rotation (corresponding to a shear- 
inversion in RPCF) occurs in the central region for Rq Ri 0.03 as detailed in figure 
me). This inversion may result from an efficient momentum transport by the large-scale 
circulation moving fast fluid from the inner wall outwards and slow fluid inwards. The 
maximum in the mean vortex torque G around Rq, ~ 0.03 indeed supports this picture, 
cf. figure [TTJ 

For the low-curvature TCF with ry > 0.9, the flattest w-profile (i.e. maximum in rdr (w)) 
and the narrow torque maximum occur for al most the same rotation numb er (i?n ~ 0.03 
versus 0.02) which extends an observation in lOstilla-Monico et al. ( 2014a l for rj ^ 0.909 
to larger radius ratios and to the limiting case of RPCF, where the flattest ita;-profile 
is linked to the narrow maximum in Nuu- This connection also exists for rj = 0.8, 
thus, supporting the independence of the narrow maximum at Rq = 0.02 from the 
detachment maximum at Rq = 0.05. For ry = 0.5 and ?y = 0.71 the Rn-value of the flattest 
angular velocity profile does not coincide with the location of the torque maximum, 
which is another evidence that this detachment maximum differs from the na rrow torque 
maximum observed for rj > 0.9. In contrast, lOstilla-Monico et al\ (|2014Qh found the 
flattest w-profile for rj = 0.5 at Rq = 0.25 and for rj = 0.714 at Rq = 0.12 which agree 
with the corresponding torque-maximizing rotation numbers. However, the difference 
might result from the fact that their values are extrapolated from simulations at Res < 
5 X 10^ lower than the Res = 2 x 10"^ used here. 

The gradients of the w-profile (ua;-profile in RPCF) are relevant when analysing the 
momentum transport in the centre and the torque (force) needed to drive the walls. 
They measure the diffusive part of the current Jc (T“, respectively) which is mediated 
by viscosity, i.e. the last term in equations (12.121) and (12.151) . A flat profile means that the 
convective current Jc,c = {urC) {Jx,c = (uyUx) in RPCF) accomplishes all momentum 
transport in the middle. 


5.3. Turbulent fluctuations 

Finally, we study the effect of the mean system rotation {Rq) on the characteristics of the 
convective momentum transport by turbulent fluctuations. We demonstrate that not all 
deviations from the laminar state contribute to the net convective momentum transport 
but only the turbulent fluctuations around the mean profile. This becomes apparent 
when introducing the velocity deviation from the mean profile u' — u — (it) by using the 
original average (•••) = (•• v differs from u" — u — {u)^ which was 

introdu ced in (14.1.1lto descr i be tur bulent fluctuations around the mean vortical motion, 
see also lBilson &: Bremhorsd ( 2007 ). As a result of incompressibility (12.61) the mean radial 
flow vanishes, i.e. (ur) = 0, so that = Ur- Moreover, by substituting the azimuthal 
velocity decomposition = {ufl) + into the convective angular momentum current. 


Jc,c = r {ur ufl) = r (ur) {ufl) + r {ur u'^) = r (m). u'^) , 


(5.4) 


we see that only the correlated fluctu ations n(. and u',„ co ntribute to the net convective 
transport. Inspired by the analysis of Burin et al\ ( 2010l) . we measure their amplitude 


by area-time averaged root-mean-square values and that enable the 

definition of a cross-correlation coefficient 

« o 


Crip — 


(^r-)rms 


SO that Jc,c — f (^r-)rms (^¥>)rms 


(5.5) 
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Figure 14. Fluctuation amplitude of the wall-normal velocity RPCF) in 

(a) and of the streamwise velocity ((^a;)rms) their cross-correlation coefficient 

Crip {Cyx) in (c) as a function of the rotation number Rn for Res = 2.0 x lO'^ and various radius 
ratios g. All quantities are radially average in the central region (r — ri)/d G [1/3, 2/3] where 
momentum is mainly transported by the convective current Jc,c = r (u'ru'fl) {Jx,c = {u'yu'flj). 


The corresponding quantities in RPCF are 

Jx,c = {'^y)rms (^’S) 

Consequently, we may consider increases in the convective (angular) momentum trans¬ 
port as being due to violent fluctuations of the wall-normal or streamwise velocities, or 
due to a strong cross-correlation between them. All three effects are analysed as functions 
of Rn for Res = 2.0 x 10^ in figure [T4l 

Generally, we observe for rj ^ 0.8 and RPCF that the fluctuation amplitudes and 
correlation coefficients collapse well in the unstable range with Rq > 0, which again 
indicates that low-curvature TCF and RPCF show the same rotation dependence. The 
wall-normal flow grows in amplitude when rotation sets in (i?o > 0) and is suppressed 
by strong rotation forming a maximum around Rn = 0.25 for rj ^ 0.71; this nearly 
coincides with the broad torque maximum at Rq = 0.2 for rj ^ 0.8, and it resembles the 
maximum observed in the torque G 2 D from streamwise invariant simulations (cf. figure 


Cyx - 


K <} 
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[n]). The co-occurrence of the maximum wall-normal flow with the broad peak supports 
the picture that the broad torque maximum is associated with strong vortical motion. 
While the co-occurrence holds for the plateau next to the actual maximum for 77 = 0.71 
(figure [TTb l . the correlation between the maxima is less clear for rj = 0.5. 

On the other hand, the streamwise velocity fluctuations are strongest without system 
rotation, forming a sharp maximum around Rq — 0, except for the high-curvature TCF 
with rj = 0.5 (figure ITlbl. Even though, a corresponding peak in the torques is missing 
since the narrow torque maximum at Rq = 0.02 clearly drops off towards Rq = 0 (cf. 

figure inb)- 

However, for rj ^ 0.9 the maximum in the cross-correlation coefficients at Rq ~ 0.03 
(figure I14h l comes closer to the narrow torque maximum suggesting that the latter is 
associated to a very efficient (angular) momentum transport by moderately strong but 
correlated vortices. Moreover, the peak in the cross-correlation coefficients coincides with 
the maximum in the torque G due to the streamwise-averaged mean flow (cf. figure fTTl) . 
The maximum in Crtp at Rq = 0.15 and at Rq = 0.05 for ry = 0.5 and rj = 0.71, 
respectively, does not coincide with the torque maximum for the corresponding radius 
ratio indicating that the latter differs from the narrow maximum for rj > 0.9. 


6. Summary and conclusions 

By describing TCF in the rotating reference-frame as proposed by iDubrulle et al. 
( 2 nn, 5 h . we were able to study its limit to RPCF as the radius ratio 77 approaches 1 . 
In this framework, shear and mean system rotation are quantified by the same parame¬ 
ters Res and Rq in both systems and curvature effects measured by Rc disappear with 
77 —>■ 1. The data collapse observed for 77 ^ 0.9 shows that Res and Rq are the appro¬ 
priate parameters in which to describe the transition from TCF to RPCF. Moreover, we 
analysed the analogy between the angular momentum transport in TCF and the linear 
momentum transport in RPCF by utilising the correspondence between the currents Jc 
and which are connected to the torque G and the force F in the respective geometries. 

As a consequence of curvature, the TCF for counter-rotating cylinders features a radial 
partitioning of stability that results in enhanced intermittent fluctuations in the outer 
partition. For a constant shear Res and changing system rotation, the occurrence of the 
intermittency has been linked to the torque maximum that forms as a f unction of Rq for 
77 = 0.5 and 0.71 ( van Cils et al. 2012 : lBrauckmann fc EckhardtJl2013(J l. Our simulations 
show that the onset of enhanced outer fluctuations agrees in the range 0.5 ^ 77 ^ 0.9 
wi th the prediction Ur, for the occu rrence of intermittency in the outer region derived 
by IBrauckmann fc Eckhardt] (I2013M) . Furthermore, we demonstrated that the stability 
difference between the inner and outer partition, and therefore the fluctuation asymmetry, 
decrease with declining curvature and eventually disappear as 77 1. As a consequence, 

also the detachment torque maximum vanishes for 77 > 0.9. 

For this low-curvature TCF, we identified two new torque maxima, a broad one at 
Rq = 0.2 that exists for all studied Res and a narrow one at Rq = 0.02 that emerges as 
Res increases to 2 x 10'^. The broad torque maximum is accompanied by a strong vortical 
flow as evidenced by the large wall-normal velocity and can be reproduced in streamwise 
invariant simulations which also points to streamwise vortices causing a high angular 
momentum transport. On the other hand, the narrow torque maximum coincides with 
an effective mean flow where the angular momentum transport is maximized not due to 
large-amplitude vortices but because of correlated vortices resulting in an effective con¬ 
vective transport. The momentum transport in RPCF features the same characteristics 
and also shows the two maxima. 
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Moreover, our simulations at Res = 2 x 10"^ reveal that Nusselt numbers, profile 
slopes, and fluctuation amplitudes collapse well for t] > 0.9 when plotted as a function 
of i?n and agree with the corresponding quantities in RPCF. These results empirically 
demonstrate the convergence of TCF with small but finite curvature to RPCF. In this 
limit, the dependence on the system rotation parametrized by Rq becomes universal in 
both systems. In addition, the angular momentum profiles in TCF show this universal 
behaviour since they collapse remarkably well for different radius ratios rj ^ 0.5 as long 
as no stabilised outer region occurs. Furthermore, we demonstrated that for the unstable 
flow {Rq > 0) the turbulent mean profiles in TCF and RPCF feature gradients in the 
centre that conform with neutral stability. Rayleigh’s inviscid neutral stability criterion 
has an analogous counterpart in RPCF and implies that the central profile slopes rdr (w) 
and analogously —dy{ux) approximate —Rq. In the low-curvature TCF and in RPCF, 
the flattest w-profile (or Ux-profile) and the narrow torq ue maximum appear a t almos t 
the same i?n value as also found for smaller t] ^ 0.909 bv lOstilla-Monico et al. (2014a). 
However for ry = 0.909, these authors observed only one maximum in the torque for 
Res 10® to 10® which appears at = 0.04 close to our narrow maximum at Rn = 0.02 
and shows the same signature in the profile gradient. This fact and the stronger growth 
of the narrow torque maximum with Res suggest that it will outperform the broad 
maximum as the shear increases further. Further simulations or experiments for low- 
curvature flows with T] > 0.9 and Res beyond 2 x 10^ could provide welcome tests of this 
prediction. 

The study of the transition from TCF to RPCF has turned out to bring more changes 
in the flow structures and the turbulent characteristics than anticipated. Perhaps the 
most intriguing observation is the transition from a single maximum in torque to two 
maxima for r] > 0.9. We could here characterize the two maxima by their relation to 
prominent features of the turbulent flow and momentum transport. While this provides 
some insights into their origin, a better physical understanding of the maxima and an 
explanation for their remarkably universal position in Rq would be desirable. 
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